Theory of Disordered Itinerant Ferromagnets II: Metal-Insulator Transition 
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The theory for disordered itinerant ferromagnets developed in a previous paper is used to construct 
a simple effective field theory that is capable of describing the quantum phase transition from a 
ferromagnetic metal to a ferromagnetic insulator, ft is shown that this transition is in the same 
universality class as the one from a paramagnetic metal to a paramagnetic insulator in the presence 
of an external magnetic field, and that strong corrections to scaling exist in this universality class. 
The experimental consequences of these results are discussed. 
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I. INTRODUCTION 



It is well known that interacting electrons in the pres- 
ence of quenched disorder at zero temperature form a 
disordered Fermi-liquid or paramagnetic metal state that 
shows, with increasing disorder, an instability against the 
formation of an insulator. This Anderson-Mott transi- 
tion (PM to PI in the schematic phase diagram shown 
in Fig. is believed to be the metal-insulator transi- 
tion observed in doped semiconductors and other disor- 
dered electron systems, ancL-i|t| has been studied theoret- 
ically in considerable detaiLuu Similarly, with increasing 
exchange interaction, the Fermi liquid state is unstable 
against the formation of long-range ferromagnetic order 
(PM to FM in Fig. |l|). This quantum phase transition 
has also been studied, both with and without quenched 
impurities. □ In the Fermi-liquid phase, the PM-PI transi- 
tion is preceded by nonanalyticities of various observables 
(e.g. the conductivity, the tunneling density of states, 
the spin susceptibility, etc.) as functions of wavenum- 
ber, frequency, or temperature. These nonanalyticities 
are often referred to as "weak-localization effects" . They 
are caused by soft modes, viz. diffusive particle-hole ex- 
citations ( "jdiffusons" ) , and can be studied in perturba- 
tion theoryB The diffusons are known to drive the metal- 
insulator transition, at least near two-dimensions, which 
is the lower critical dimensionality for this transition. 
The analogous soft-mode effects in the metallic ferromag- 
netic state have recently been investigated in Ref. ^| (to 
be referred to as (I)), but the quantum phase transition 
that must occur from a ferromagnetic metal to a ferro- 
magnetic insulator upon increasing the disorder (FM to 
FI in Fig. g) has never been considered. 

In this paper we address the latter problem. In par- 
ticular, we derive and analyze an effective field theory 
that is capable of describing the disorder driven transi- 
tion from a ferromagnetic metal to a ferromagnetic in- 
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FIG. 1. Schematic phase diagram of disordered, interacting 
electrons, showing paramagnetic metal (PM), paramagnetic 
insulator (PI), ferromagnetic metal (FM), and ferromagnetic 
insulator (FI) phases. M denotes a multicritical point. 



sulator. An interesting theoretical question that arises 
in this context is the role of the Goldstone modes that 
occur due to the broken spin rotational symmetry, i.e. 
the spin waves. Since they constitute soft modes in ad- 
dition to the diffusons, one would a priori expect them 
to influence the critical behavior. It was shown in (I) 
that the Goldstone modes, while they contribute to the 
leading frequency nonanalyticity of 0(fl( d ~ 2 ^ 2 ) in the 
conductivity, yield a prefactor that is of O(l), while the 
diffusons contribute a prefactor that is of 0(l/(d — 2)) 
and thus diverges as d — > 2. Since it is known that this 
singularity drives the transition near two-dimensionsJ3 it 
follows that the Goldstone mode do not contribute to the 



1 



asymptotic critical behavior. Largely as a consequence of 
this, the ferromagnetic-metal-to-ferromagnetic-insulator 
transition turns out to be in the same universality class 
as the one from a paramagnetic metal to a paramagnetic 
insulator in the presence of an external magnetic field. 
We had strong corrections to scaling for this universality 
class. □ 

Another important motivation for the present study 
is the recently observed apparent metal-insulator tran- 
sition in Si MOSFETs and other two-dimensional (2-d) 
electron systemsp which contradicts the orthodox theo- 
retical results that predict an insulating state in d = 2 
even for arbitrarily weak disorder .EMa Since it is known 
that magnetic fluctuations have a tendency to increase 
the conductivity in or close to two-dimcnsions,EO it is 
conceivable that there might be a ferromagnetic metal- 
lic phase at small but nonzero disorder in d = 2. We 
find that this is not the case, which rules out a possible 
mechanism for a metal- insulator transition in d — 2. 

This paper is organized as follows. In Sec. [0] we briefly 
recall the general Q-matrix field theory for itinerant fer- 
romagnets that was developed in (I). On the basis of this, 
and using the perturbative results of (I), we construct an 
effective theory for the most relevant soft modes in the 
system. This theory takes the form of a generalized non- 
linear a model. In Sec. Ill we show that this model de- 



scribes a ferromagnetic-metal-to-ferromagnetic-insulator 
transition in d > 2, and we calculate the critical behavior 
at this transition in a d = 2 + e expansion. In Sec. |iy| we 
conclude with a discussion of our results. 



II. EFFECTIVE FIELD THEORY FOR 
DISORDERED ITINERANT FERROMAGNETS 



(2.1c) 



r,i— 



and analogously for A. Here tq = sq = I2 is the 2x2 unit 



matrix, and tj = ~Sj — —iaj, (j = 1,2,3), with CTi.2,3 
the Pauli matrices. In this basis, i = and i — 1,2,3 de- 
scribe the spin singlet and the spin triplet, respectively. 
An explicit calculation reveals that r = 0, 3 corresponds 
to the particle-hole channel while r = 1 ,2 describes the 
particle-particle channel. In Eq. ( 2.1a ), Tr denotes a 



trace over all degrees of freedom, including the continu- 
ous position variable, while tr is a trace over all those 
discrete indices that are not explicitly shown. For the 
disorder part of the action one findsEj 



A dis [Q] 



1 



dx tr (Q(x))' 



(2.2) 



irNpT , 

with r the single-particle scattering or relaxation time.B 
The electron-electron interaction A- m t is conveniently de- 
composed into four pieces that describe the interaction 
in the particle-hole a&d particle-particle spin-singlet and 
spin-triplet channelsE3 We will need only the particle- 
hole channel, and thus write 



A int [Q] = A^ +A- 



1 W 



(2.3a) 



A' 



(«) 



7T(s) 



dx ^ (-1)' 

r=0,3 



E E 



ni,Ti2,m a 



x [tr ((r r ® So )QZn 1+m ^))} 



x [tr ((r r ® s ) Q 



aa 

n2+m,7i2 



W)] 



(2.3b) 



A. Q-matrix theory 

In (I) it was shown that disordered itinerant ferromag- 
nets are described by the following action, 

A[Q, A] = A dis + Ant + \ Tr In (g o 1 - iX) 

+ J dx tr (A(x) Q(x)) . (2.1a) 

Here 

Gq 1 = —d T + c?x/2m c + n , (2.1b) 

is the inverse free electron Green operator, with d T and d x 
derivatives with respect to imaginary time and position, 
respectively, m c is the^ electron mass, and \i the chemi- 
cal potential. Q and A are matrix fields that carry two 
Matsubara frequency indices n and m, and two replica 
indices a and (3. The matrix elements Q"m and A"^ are 
spin-quaternion valued. They are conveniently expanded 
in a basis 



^~/^x £(-ir E EE 

r— 0,3 ni,ri2,m a i—1 

x [tr ((T r ® Sl )QZ ni+m (x))} 

x [tr ((r r (g) Si )Q^ +mi „ a (x))] . (2.3c) 

Here r( s ) > and r' 4 ' > are the spin-singlet and 
spin-triplet interaction amplitudes, respectively. T^' is 
responsible for producing magnetism. 

As was shown in (I) , the Goldstone modes do not con- 
tribute leading singular terms to the conductivity, or to 
any of the two-point vertices, in perturbation theory up 
to one-loop order as d — > 2 in the metallic ferromagnetic 
phase. It was further shown that they do not contribute 
to the renormalization of the DOS, and hence they do 
not contribute to the wave function renormalization. It 
follows that both the heat and the charge diffusion con- 
stants do not carry any singular renormalizations due to 
Goldstone modes. If we assume that this signals the ab- 
sence of Goldstone mode effects on the critical properties 
near the metal-insulator transition as well, we can ignore 
the Goldstone modes for the purpose of constructing an 
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effective field theory for the soft modes that drive the 
metal-insulator transitionO Furthermore, we can ignore 
the particle-particle channel, which is massive in a sys- 
tem with a nonvanishing magnetization. Accordingly, we 
drop both the particle-particle channel (r = 1, 2) and the 
transverse spin-triplet channels (i = 1, 2) from our model 
definition. For the remaining soft modes, we will now con- 
struct an effective theory by generalizing the procedure 
followed in Ref. O. 



B. Soft and massive modes 

Let us briefly recall the basic philosophy behind the 
derivation of a nonlinear cr-model in Ref. [l2], which in 
turn was based on the work by Schafer and WegnerEj on 
nonintcracting electrons. First one realizes, by means of 
a Ward identity, that the soft modes are given by the 
matrix elements Q nm with nm < 0, while the Q nm with 
nm > are massive. This remains true in the present 
case except for the Goldstone modes, which we can ne- 
glect for our purposes. Next one block-diagonalizes the 
matrix Q in frequency space. Algebraic arguments show 
that the most general Q can be written as 



Q = SPS- 



(2.4a) 



Here S is a matrix that represents an element of the coset 
space USp(8An,C)/USp(4An,C) x USp(4An,C), and P 
is block-diagonal in Matsubara frequency space 



P = 



P> 






p< 



(2.4b) 



where P > and P < are matrices with elements P nm where 
n, m > and n, m < 0, respectively. It is further conve- 
nient to define a transformed field A by 

A(x) =5- 1 (x)A(x) l S(x) , (2.5) 

and to write the action in terms of these variables, 

A [S, P, A] = A dis [P] + AntiSPS- 1 } 

+— Tr In (Gn 1 — iSAS^ 1 ) 



+ J dxtr (A(x)P(x)) 



(2.6) 



The next step is to expand S, P, and A about their 
respective saddle-point values, which we d enote by (S), 
(P), and (A), respectively. From Sec. II B in (I) we have 



(S) = 1 <g> r , (2.7a) 



P)u = 5 12 ^ t ( T ° s °) &i (P) 
p 

+ (r 3 ®s 3 ) ?m(p)] , (2.7b) 



S12 (t <8> so) 



ttA f t V 



77 G ^ (P) 



S12 (r ®.s ) 2iT^tJ2^° I^£ m (p) 



+ $12 (t 3 ® s 3 ) 



itN f t V 



S12 (r 3 ® S3 ) 2 l r( t )T^ e ^° I^ m (p) , 



(2.7c) 



with Q and J- from (I) Eqs. (2.13). In the popular approx- 
imation that replaces the wavevector sum_over a Green 
function by an integral over £ p = p 2 /277i e ,E3 we have 



1 \ -s . . — ITT 

— Snip) ~ — N F sgnw n 



— ITT 
2~ 



t^E^«(p) 







(2.8a) 
(2.8b) 



In this approximation,0 we can write Eq. ( [2.7b] ) as 

(P)i2~iri2 , (2.9a) 

with 

7Ti2 = S 12 (t <S> s ) sgnw ni , (2.9b) 

with u) n = 2irT(n + 1/2) a fermionic Matsubara fre- 
quency. For our purposes this approximation will be suf- 
ficient for reasons that were explained in detail in Ref. 

We now write 



P = {P) + AP , A = (A) + AA 



(2.10) 



and expand in powers of AP, AA, and deri vativ es of S. 
Let us first consider the Tr In term in Eq. (2.6). Using 
the cyclic property of the trace, we can write it in the 
form 

Tr HGo 1 -SiAS- 1 ) = Tr h^S" 1 G^S - iA) 
= Tr ln(G sp ) _1 ) + Tr In [ 1 + G sp S^^S) 

+— G sp S- 1 (V5) V + G sp S-\V 2 S) 
m 2m 



-G sp i(AA) 



with 



G s 



(Go 1 



i<A»- 



(2.11a) 



(2.11b) 



the saddle-point Green function. This is formally the. 
same expression as in the absence of ferromagnetism,E3 
only the saddle-point Green function is more compli- 
cated. In particular, the transformation matrix S ap- 
pears only in conjunction with some derivative and is 
therefore soft, while the fluctuations AA are massive. 
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Expanding the second term on the right-hand side, the 
simplest contribution is the one involving the time deriva- 
tive, 

Tr [GspS^idrS)] 

= Jdxti [ifi<S(x)G sp (x = OjS-^x)] 



t:N f 



dx tr [Q Q(x)] + 0(n 2 Q) . (2.12) 



Here 



12 



(r ® S ) <5l2 fit!! 



(2.13a) 



is a frequency matrix with £l n — 2nTn a bosonic Mat- 
subara frequency, and 

Q(x) = S(x)wS- 1 {x) , (2.13b) 



with p-.hom Eq. (2.9b). Here we have made use of Eqs. 
l ^ This is the same result as the one obtained in the 



absence of ferromagnetism.E3 

We now turn to the gradient terms. It is convenient to 
define a matrix valued <i-dimcnsional vector field 



s(x)=5- 1 (x) (V5)(x) 



(2.14) 



and to expand in powers of s. The term linear in s van- 
ishes for symmetry reasons. To 0(s 2 ), bo th the next-to- 
last term on the right-hand side of Eq. ( [2. 11a ) and the 
square of the preceding term contribute. So far our gra- 
dient expansion has been completely general. In order 
to evaluate the terms of 0(s 2 ), we now remember that 
we can neglect the spin waves for the purpose of deriv- 
ing a soft-mode transport theory, i.e. we have dropped 
the channels r = 0, 3, i = 1, 2 in the spin quaternion ex- 



pansion, Eq. (2.1c). Furthermore, it is well known that 
the Cooper channel (r = 1,2) is massive in an external 
magnetic fieldjja and the same is true in a ferromagnetic 
state. Consequently, the only spin-quaternion degrees of 
freedom present in S12 are to,3 and So,3, and S12 com- 
mutes with T3 ® S3. This simplifies the evaluation of the 
gradient terms substantially, and we obtain 

Tr G sp S- 1 (V 2 S) - - Tr (G sp S' 1 (VS) V) 2 = 



EE 

12 q 



dx 



»7fa,«(q)fr (*i2(q)*ai(-qC 



+Vi2,iM) tr (( r 3 ® s a) s i2(q) 4i(-q))] . (2.15a) 



where 

rf = (r,+ + r,-)/2 
with 

,± 



1 



+— %5j^ 1 (q)^ 2 (q) 



(2.15b) 



(2.15c) 



with the Green functions defined in (I) Eq. (2.13c). 
Equations ( 2.15| ) are generalizations of the corresponding 
expressions in the absence of ferromagnetism.113 



C. The nonlinear a model 

The remaining steps in the derivation of an effective 
field theory proceed in analogy to Ref. [l^. In particular, 
we integrate out the massive modes P and A in tree ap- 
proximation, i.e. we neglect all fluctuations AP and AA. 
Tj^~ can be related to the conductivity in self-consistent 
Born approximation of a system whose chemical potential 
has been shifted from its value for nonmagnetic electrons 
by ±A = ±T^'M/ fJ,B- Here M is the magnetization in 
Stoner approximation, and /is is the Bohr magneton (see 
(I) Eq. (2.15)). Denoting these conductivities by <Jq, and 
defining the bare coupling constants 



1/G - - m (a+ + a ) 

I/G3 = ^m(4 -fo") 
H = nN F /8 , 



we obtain for the effective action 



(2.16a) 

(2.16b) 
(2.16c) 



A=—Jdxtr [VQ(x)] 2 + 2H J dx tr [fi Q(x)] 
J dx tr ((r 3 ® s 3 )[VQ(x)] 2 ) + A int [Q] 



(2.17) 



Here Q — Q — ir with ir t he matrix defined in Eq. ( 2.9b| ) 
and Aint from Eqs. (2.3). We recognize this action as 
the generalized nonlinear a model for disordered inter- 
acting electrons^ augmented by the term with coupling 
constant I/G3 that is proportional to the magnetization. 

It turns out that the bare action, Eq. (2.17), is not 
sufficient to completely describe the effects of magnetic 
long-range order, even if one ignores the spin waves as we 
did. As we will see, under renormalization two additional 
terms are generated. One is a frequency coupl ing t hat is 
analogous to the second gradient term in Eq. ( [2.171 ) 1 an d 
the other is an electron-electron interaction term that is 
not present in nonmagnetic systems. We therefore need 
to add these terms to our action. Denoting the respective 
coupling constants by H 3 and K 3 , we obtain our final 
result for the effective action, 

A = ^ f dx tr [VQ(x)] 2 + 2H j dx tr [fi Q(x)] 



^-Jdxtr ((r 3 ®s 3 )[VQ(x) 



2G 



2H 3 dxtr (T 3 ®s 3 )nQ{x) + A hlt [Q} . (2.18a) 



Here ^4i nt is given by Eqs. (2.3) plus the extra term. In- 
troducing new interaction amplitudes K s = — 27rr( s \-&nd 
K t — 27iTW to comform with notation used earlier, a we 
write 



Ant [Q] = 4nl [Q] + [Q] + 4% IQ] 



(t). 



i(3)r 



(2.18b) 
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-ttT 



K s / dx 2_j ^qiq 2 3a ia3 #1-2,4-3 
J 1234 



x ^H r tr [(r r ®«o)Qi2(x)] 

r 

xtr [(tv®s )Q 34 (x)] , (2.18c) 

Ant[<9] = J lf '/ (ix I! ^i"2 5 aia3 $1-2,4-3 
1234 

X HH rtr [( T r®S3)Ql2(x)] 

r 

xtr [(T r <g)s 3 )Q 34 (x)] , (2.18d) 



AntlQ] = -47rTi^ 3 / dx ^ S <* ia2 S aia3 Si- 

-EE 



2,4-3 



1234 

m rs ,ij ~ r Q 12 (-x)lQ 34 (x) , (2.18e) 



where 



itr(r 3 T r rj') tr (s 3 s 4 s}) 



(2.18f) 



(This is the matrix that was denoted by m 03 in (I).) 



Finally, we note that Q as defined in Eq. ( 2.13b ) obeys 
g 2 (x) = l , Q+ = Q , trQ(x)=0 . (2.19) 



Equations ( 2.18] . 2~T^ ) represent the analog for itinerant 
ferromagnets of the nonlinear a modelo for paramagnetic 
electron systems. 



D. The metallic fixed point 

Before we proceed to use the a model to study the 
quantum phase transition from a ferromagnetic metal to 
a ferromagnetic insulator, let us ascertain that the model, 
with some correction terms, actually describes a metallic 
ferromagnetic phase in some parts of parameter space. 
Since we will want to approach the transition from this 
phase, its existence within the model is obviously a nec- 
essary condition for our program to be viable. 

This task is very simple, since it proceeds in exact anal- 
ogy to the demonstration in Ref. [L2| that the model with 
I/G3 = H3 = K3 = has a stable Fermi-liquid fixed 
point. This is because the power counting procedure used 
to prove the existence of a stable fixed point does not de- 
pend on structural details like the presence of extra r 
and s matrices in the various terms of the action, while 
such details are the only difference between the current 
model and the one considered in Ref. [l^. Accordingly, we 
parameterize Q in terms of a matrix q with elements q nm 
whose frequency labels are restricted to n > 0, m < 0, 



Q 



93 



q 



and expand S in powers of q, 



5 = 1 



1 ( 



' T0+ 2U t 







0(q 2 ) 



(2.21) 



As in Ref. [lj, we assign scale dimensions to <?(x), 

[«(x)]=(d-2)/2 , (2.22a) 
and to the fluctuations of the fields P and A, 



[AP(x)] = [AA(x)] = d/2 



(2.22b) 



Here the scale dimensions [. . .] are defined such that the 
scale dimension of a length L is [L] = — 1. The fixed point 
action then consists of the nonlinear a model action, Eq. 
( |2.18a ), expanded to 0(q 2 ), plus the corrections bilinear 
in AP and A A that arise from Eq. (pij|). All other terms 
are irrelevant by power counting. The arguments show- 

fthis are exactly the same as the ones given in Ref. 
and need not be repeated here. The correlation func- 
tions for this Gaussian action are simply related to the 
Gaussian propagators of Sec. Ill in (I). We will explicitly 
determine them in Sec. |l| below. This will show that the 
fixed point action really describes a disordered itinerant 
ferromagnet. 

In contrast to Ref. |l2], however, we cannot discuss the 
leading corrections to scaling near the stable metallic 
fixed point within our current framework. The reason is 
our having neglected the transverse spin-triplet channel 
that contains the Goldstone modes. While the latter are 
not expected to influence the leading scaling behavior at 
the critical fixed point for the reasons pointed out above, 
they do contribute to the corrections to scaling near the 
metallic fixed point, as indicated by their contribution 
to the leading nonanalytic frequency dependence of the 
conductivity that was studied in (I). 



III. METAL-INSULATOR TRANSITION ON THE 
BACKGROUND OF FERROMAGNETISM 

In this section we perform a o ne-lo op renormalization 
of the nonlinear a model, Eqs. (2.18). We first do this 
for general parameter values, which leads to rather com- 
plicated flow equations. They contain a fixed point that 
corresponds to the known critical fixed point fojr nonmag- 
netic electrons in an external magnetic field.Ej We then 
linearize about this fixed point and show that it is per- 
turbatively stable with respect to the additional terms 
in the action that represent the presence of a nonzero 
magnetization. 



A. Parametrization, and Gaussian order 



In order to set up a loop expansion we use the param- 



(2.20) eterization for the matrix Q that is given by Eq. (2.20) 
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Note that thi s par ameterization builds in the constraints 
given in Eq. ( 2.19 ). It is the matrix analog of the usual 
(it, 7?) parameterization of the O(N) vector nonlinear a 
model.Eij The loop expansion now proceeds as an expan- 
sion in powers of q. To Gaussian order, we obtain 



Ai0) = l^EE^12(P)^ M 12,34(P)i934(-P) 



V 



p 1234 



where the Gaussian vertex is given by 
%M 12M (p)= 6 13 5 2i %M^(p) 



+<5l-2,3-4 5 aia2 S aia3 2~kT K r 



(3.1a) 



(3.1b) 



with 



%M l S(p) = S ra 5ij - (p 2 + Gim ni _„ 2 ) 



+m r 



G, 



p 2 + G 3 H 3 n ni _ n2 ) 



(3.1c) 



and 



K rStij = S rs S i: j (S M K S + §izK t ) + m rs ,ij K 3 , (3. Id) 



with m rs ,ij from Eq. (2.181). 

The Gaussian propagator can be determined by the 
same methods that were employed in Sec. Ill of (I). We 
find 

C r g 12 (k)^ 34 (p))(°) = i5(k + p);iM- 2 1 > 3 4 (p) , (3.2a) 
where M _1 has the structure 

ra M 12,34(P) = ^3 ^34 [S rs $ij ^m-n 2 (p) 

+01-2, 3—4 ^a\a.2 ^a±a^ \prs $ij C ni —ri2 (p) 

+m rs ,ij -Dni-na(p)] • (3.2b) 

To specify the propagators A, B, C° = C s , C 1 ' 2 ' 3 = C\ 
and D, wc define 



a = a„(p) - (p 2 + GHVt n )/G , 
b = b n (p) = (p 2 + G 3 H 3 Q n )/G 3 , 



and 



N = N n {p) = (a 2 - b 2 ) [a 2 - b 2 - 2bK 3 Q n 
+a(K s + K t )fl n - K 2 n 2 n + K s K t n 2 n ] 

In terms of these quantities, we have 

A n (p) = a/(a 2 -b 2 ) , 
B n (p) = -b/(a 2 + b 2 ) , 
-2vrT 



(3.3a) 
(3.3b) 



(3.3c) 



(3.4a) 
(3.4b) 



C'(P) 



N 



[a 2 K s + b 2 K t + a(-2bK 3 - K$il r . 



-K s K t n n )] , 



-2ttT 



C* (p) = — ^- [a 2 K t + b 2 K s + a(-2bK 3 - K 2 VL n 



-K s K t Q n )] , 



2nT 



D n (p) = =j=- [~a 2 K 3 + ab{K s + K t ) - b{bK 3 + K 2 VL n 
-K„KSUS\ ■ 



(3.4c) 
(3.4d) 

^3 

(3.4e) 



B. Perturbation theory to one-loop order 

We now proceed to perform a one-loop renormaliza- 
tion of the theory. We do this by renormalizing the two- 
point vertex ^ S M 12 34 , Eqs. (3.1). This procedures proves 
the renormalizability of the theory to one-loop order, i.e. 
it makes sure that no coupling constants in addition to 
the ones present in the bare theory are generated under 
renormalization. Wc also need to determine the wave- 
function renormalization. This we do by considering the 
one-point vertex function T^> = (Q)^ 1 . 



1. One-point vertex 

Let us first consider the one-point propagator (Q). 
To one-loop order, the only diagram that contributes 
is shown in Fig. |^. In spin-quaternion space, there are 




<Q>=1 + 



FIG. 2. Perturbation theory for (Q) to one-loop order. 

two nonvanishing matrix elements of (Q), viz. (qQ) and 
( 3 Q). These expectation values are diagonal in both fre- 
quency and replica space. Their inverses constitute one- 

I3 (fi n ), resp ective ly. A simple calculation using the 



point vertex functions that we denote by T (f2 n ) and 

Tg (£!„), resp ective ly. A 
results of Sec. Ill A yields 



Here we have defined the integrals 

p l—n 
1 00 



(3.5a) 
(3.5b) 

(3.6a) 
(3.6b) 
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2. Two-point vertex 

We now turn to the two-point vertex I^ 2 ) , whose Gaus- 
sian approximation is given by Eqs. (3.1). To one- loop 
order, we write 

.(2) 



rs L 12,34 



(P) 



13 M 

rs lvl 12,34 ^ 



(p)+£(5M) 12 , 4 (p) . (3.7) 



There are two topologically distinct diagrammatic con- 
tributions to SM , which are shown in Fig. |^. They arise 





FIG. 3. One-loop contributions to the two-point vertex. 

from quartic terms, i.e. terms of 0(q ), and cubic terms, 
i.e. terms of 0(q 3 ), respectively, in the expansion of the 
action in powers of q. An evaluation of the diagrams 
is straightforward but very tedious, since the two topo- 
logical structures can be dressed in many ways with the 
various indices carried by the g-field. The calculation re- 
veals that the one-loop contributions can be grouped into 
three distinct classes: (1) Quartic contributions that are 
logarithmically divergent as f2„ — ► in d = 2, (2) Cubic 
contributions that have the same degree of divergence, 
and (3) contributions, both quartic and cubic, which indi- 
vidually diverge more strongly ('superdivergent terms'), 
but combine to yield again terms that are only logarith- 
mically divergent. Denoting the contributions of these 
three classes to the one-loop renormalization of by 
5M<- 4 \ 8M&\ and 5M^ sd \ respectively, we find for the 
first of these classes 

j£(M0$34(l>) = *13 ^24 UsSi~ (p 2 + GHSl N ) 



x[l['\Sl N )+I®({l N )] 



+*r.*« (p 2 + G 3 H 3 n N ) 2 i[ 3 \n N ) 



i 



+m 



8G 3 

' v ~8G ^ 2 ^^^ N ^ (P* N ) 



8G- 



( P 2 + G 3 H 3 n N ) [i[ s) (n N ) + i[ t) (n N )] 



-^l-2,3-4^a 1 Q 2 5a 1 Q 3 ^- T |(5 rs 6 4 ., [(K„ + K t ) J(Q N ) 

+2K 3 j^(n N )] 

+m rs , ij [2K 3 J(n N ) + (K s + K t ) (n N )] X , (3.8a) 



with l[ s '*' 3) given by Eqs. Q, 



j( 3 )(n JV ) = i^B JV (p) , 



(3.8b) 
(3.8c) 



and N an external frequency (e.g., N — ni — n-i). For 
the second class we obtain 

^(<5Af)^3 4 (p) = -<Sl-2,3-4<5a 1 a 2 <5Q 1 Q 3 ( 7rr ) 2< | <5, .A/ 

'2\ As) . 



x [(S M K 2 + 6 l3 K s K t + Kl) li s> (fi w ) 
+Mt + S i0 K s K t + K 2 ) jf^jv) 
+ {8 M K S K 3 + K S K 3 + 8 l3 K t K 3 ) 4 s) {Q N ) 
+ {8 lQ K s K 3 + 6 a K t K 3 + K t K 3 ) /<*> (Q N )] 



+m rSiij [(S M K 2 + 5 j3 K s K t + K 2 ) 4 s) (fl 



N) 



+{5 i3 Kl + 6 l0 K s K t + K 2 ) I { i\n N ) 
+ {8 lQ K s K 3 + K S K 3 + 5 l3 K t K 3 ) I [ f' \(l N ) 

+ (S M K S K 3 + 5 a K t K 3 + K t K 3 ) I^^n)] 

(3.9a) 

Here we have defined integrals 

ii s ' t) (n N ) = ±j2Y,[Mp) 2 +B l ( P f 

p l—N 

+A ; (p) I Cl s ' t] (p) + Bx (p) / Di (p)] , (3.9b) 

1 OO 

4 s ' t} (n N ) = -J2J2[ 2A i(p)Bi(p) 

p l — N 

+B l (p)lC^ t) (p)+A l (p)lD l (p)] , (3.9c) 
Finally, for the third class we have 

%(SM)[ s 2 %(p) = S 13 5 2i I [4 S) ( P , n N ) 

+4 t \p,n N )] 

+m rs , ij [4 s) (p,n N )+I 3 :t) (p,n N )]} , (3.10a) 
with integrals 

oo , 

4 s < 4) ( P , a N ) = - W k ) c / ( k ) + w D ^ 

k 1=1 ^ 

+2nTK s>t Mk) + 2nTK 3 B,(k) 

- (2tt T fK 2 1 1 A l+N (k + p) [At (k) + I C, (M) (k)] 

- (2tt TfK s K t I B l+N (k + p) [B t (k) + I D l (k)] 
-(2tt T) 2 if 2 i A ;+W (k + p) [At (k) + i Cf ' t} (k)] 
- (2tt T) 2 if 2 i B l+ N (k + p) [B t (k) + i A (k)] 
-(2TTT) 2 2K 3 K s , t l A l+N (k + p) [25/ (k) + Z A(k)] 



-(27rT) 2 2^ 3 ^ s , t I S /+7V (k + p) i C, (S|4) (k) 



(3.10b) 



OO 

k /=1 ^ 
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+27r T K s>t Bi (k) + 2ttTK 3 A t (k) 
- (2tt T) 2 A 2 1 Z B l+N (k + p) [At (k) + I C\ s ' t] (k)] 
-(27rT) 2 A s A t i A i+A r(k + p) [B,(k) + / D,(k)] 
-(2tt T) 2 A 3 2 i £ i+JV (k + p) [A, (k) + / C, (i,t) (k)] 
-{2-nTfKll A l+N (k + p) [B Z (k) + / D t (k)] 
-(27rT) 2 2A 3 A s , t / A ;+A r(k + p) [A,(k) + i Cf M) (k)] 

-(27rT) 2 2A 3 A s , t / B i+JV (k + p) [B,(k) + i D,(k) 

(3.10c) 

As a check, we consider the superdivergent contributions 
at zero external frequency and wavenumber. Although all 
of the individual terms are linearly divergent, an explicit 
calculation yields 



4 s) (0,0) + 4 t) (Q,0) = , 
and the same for I 3 . Therefore, 

%(5M)&%(p) = 0(p 2 ,n N ) 



(3.11a) 



(3.11b) 



with coefficients of p 2 and f2jv that are only logarithmi- 
cally divergent. The cancellation of the superdivergen- 
cies thus holds as expected (and required by, e.g., parti- 
cle number conservation and the renormalizability of the 
theory. ) 



C. Expansion to linear order in the magnetic 
coupling constants 

As is clear from the preceding subsection, the complete 
one-loop renormalization of our model is rather compli- 
cated. While it is certainly possible to determine the 
renormalization group (RG) flow equations from our per- 
turbative results, it would not be easy to analyze them 
for fixed points. At this point we therefore take a less 
general approach that is based on the following physical 
considerations. We are interested in a phase transition 
from a metallic magnetic phase to an insulating mag- 
netic phase. Physically, we expect the magnetization to 
be noncritical at such a transition. The simplest possible 
scenario is then a fixed point where the renormalized val- 
ues of the 'magnetic' coupling constants, I/G3, H 3 , and 
A3, are all zero. (More complicated possibilities we will 
come back to in Sec. |y] below.) This means that the un- 
versality class of this transition is the same as that for the 
transition from a paramagnetic metal to a paramagnetic, 
insulator in the presence of an external magnetic field.tl 
We can check this scenario by expanding to linear order 
in the three magnetic coupling constants and investigate 
the perturbative stability of the nonmagnetic fixed point. 

Accordingly, we expand the results of the previous sub- 
section, expressing the result in the form of corrections 
to the magnetic coupling constants. We use dimensional 



regularization, i.e. we perform the integrals in d = 2 + e 
to leading order in 1/e. We find for the correction to 
I/G3 to linear order in I/G3, H3, and A3, 



5(1/G 3 



G 



16G : 



f n (K s /H,K t /H) 



-|| MK./H, K t /H) + H fi 3 (K s /H, K t /H) 



(3.12a) 



Here we have defined the functions 
fn(x,y) = 9n(x) + gn(y) - 2 Ll 



x-y 



x-y 



[ftii(x) - hu(y)} 



fi2(x, y) = g 12 (x) + ffi2 (y) + — — — 

x-y 



fu(x,y) 


X 


in terms of 


9n(x) = 


-6 

X 


hn(x) = 


-2 

X 2 


012 (a) = 


1 

X 


hi 2 (x) = 


1 

x z 


and L x = 


ln(l 



xy 
x — y 



[h 12 (x) - h 12 (y)} 



x — y \x 



6 / 1 

x \ x 



x* x x* V x 
3 



X 



+ 2 In(l + ar) 



1 ln(l + x) 



(3.12b) 

(3.12c) 
(3.12d) 

(3.12e) 
(3.12f) 
(3.12g) 
(3.12h) 



Similarly, the correction to H3 is 



^-MK./H.KJH) 



-H 3 f 22 (K s /H, K t /H) + A 3 fe{K./H, K t /H) . 

(3.13a) 



with 



f2i{x,y) = 021 (x) +521(2/) - 1 1 X ~ V ( L * - L v) 



xy 



x-y 

f-n{x,y) = g 22 {x) + g 22 (y) 
xy 



x-y 
h 21 {x) - h 21 (y)} 

L X — Ly 



x-y 



h 21 (x) - h 21 (y)} 

x-y 

f23(x,y) = g 23 (x) + g 23 (y) 1 



x-y 



+ [h 21 (x) - h 21 (y)} 

x-y 



(3.13b) 



(3.13c) 



(3.13d) 
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in terms of 
321 0) : 
h 21 (x) ■■ 
922{x) - 

923(x) ■■ 



1 + E-fI + l) inci +,) 



1 



ln(l + x) 
1 



-1+ (- + 1 | InCI + .<■) 



-1 1 



ln(l + a;) 



2 x 

Finally, for the correction to K 3 we obtain 



(3.13e) 
(3.13f) 
(3.13g) 
(3.13h) 



SK 3 = ? 



^-f 31 (K s /H,K t /H) 



+#3 f 3 2(K s /H, K t /H) + K 3 f 33 (K s /H, K t /H) 

(3.14a) 



with 

^ + y 

/31 (a:, y) = x + y H [h 31 (x) 

x-y 

y 

h2{x,y) = g 32 (x) +532(2/) + -L a 

x 



h3i(y)] 

X T 

y 



(3.14b) 



[x + y? 

x-y 



[h 32 {x) - h 32 (y)} , (3.14c) 



hz{x : y) = g 33 (x) + g 33 (y) + - L x 

x 



in terms of 



h 31 (x) 
932{x) 

h 32 (x) 



+ {x + yf [h32{x) _ hM] 
x-y 



-2(1 + a;) ln(l + x) 
-2x + \n(l + x) , 



ln(l + x) 



(3.14d) 

(3.14e) 
(3.14f) 

(3+4g) 
(3.14h) 



333(2:) = l + 31n(l + a;) . 

An inspection of the integrals in Sec. Ill P3| further 
shows that all corrections to the remaining coupling con- 
stants, G, H, and K a t, to the extent that they depend on 
the magnetic coupling constants, are at least quadratic 
in the latter and hence can be neglected for our purposes, 
The 'nonmagnetic' one-loop corrections are well known,EI 
and we do not write them down again. 

We also note that I/G3 7^ is sufficient to generate 
nonzero values of H 3 and K 3 in perturbation theory, even 
if these coupling constants were not present in the bare 
actio n. Th is is the reason why we have included them in 
Eq. fl2.18aP. 



D. Renormalization group flow equations 

We now perform a RG analysis of our perturbation 
theory. We define renormalizcd coupling constants 173, 
h 3 , and k 3 by 



G 3 = K 
H 3 = Zh 3 h 3 
K 3 = Z k3 k 3 



,93 



(3.15a) 
(3.15b) 
(3.15c) 



where the Z are renormalization_cpnstants, and k is the 

arbitrary RG momentum scale.Ej We further define a 

(2) 

two-point vertex function T 3 as the 'magnetic piece' of 



the general vertex defined in Eq. (3.7), i.e. the parts 
t hat are proport ional to I/G3, H 3 , and K 3 . From Eqs. 
(|+]J33,UH3T^) we have 



r 3 2) ( P ,fi) 



G 3 



S(l/G 3 ) p 2 + (H 3 + SH 3 ) n 



+ (k 3 + 6K 3 ) n 



(3.16) 



The renormalization constants can then be determined 
from the renormalization statement 



rg^p, n ; g 3 ,h 3> k 3 -,K) = zrn P) n; g 3 ,h 3 ,k 3 ) , 



(2)/ 



(3.17) 



(2) ) 
where T 3 L is the renormalized counterpart of T]^' , and 

Z is the wavefunction renormalization. In our notation, 
we suppress the dependence of the vertex functions on 
the remaining coupling constants, G, H, K Si t, and their 
renormalized counterparts. 

It is a priori not clear that a single wavefunction renor- 
malization constant will suffice. Indeed, the existence of 



(2) 



two distinct one-point vertex functions, Eqs. (3.5), one 
being related to the density of states and the other to the 
magnetization, might suggest that one needs at least two. 
However, as mentioned in Sec. Ill C above, we do not ex- 



pect the magnetization to display leading critical behav- 
ior at the phase transition we are interested in, despite 
the fact that the magnetization has nonanalytic contri- 
butions in perturbation theory. We therefore expect the 
only wavefunction renormalization to be the one related 



the vertex T 



(i) 
> 



if^fi; 33, h 3 , k 3 ; K ) = Z if > (fl; G 3 , H 3 , K 3 



(3.18) 



To linear order in 1/33, h 3 , and k 3 , Z is given by the 
wavefunction rcnormalizatiorufor nonmagnetic electrons 
in an external magnetic field,EI 



4e 



(ls+h) 



(3.19) 



where l s> t = ln(l + 7 s> t), with j s j = k s ,t/h the renormal- 
ized counterparts of K s t /H. Our perturbative calcula- 

(2) 

tion of T 3 is then sufficient to determine the remaining 
renormalization constants. Using minimal subtraction, 
we find 



Z„ 



9_ 

8e 



/n( 7s , It) - 2(1. + h) + 2 ^ fu (is , It) 



gh 



+4^^ /l3(7a,7i) 

gh 



(3.20a) 







9_ 

4c 



Z k , = 1 + ^ 



h + lt~ j '22(7s, 7*) 

-7^-/23(78,7*) 

«3 



53^3 



/2i(7s,7t) 



(3.20b) 



2(h + h) - /33(7s,7t) - 7-/32(7s,7t) 

fc 3 



53^3 



/i3(7s,7t) 



(3.20c) 



From Eqs. (3.2C) and (3.15) it is now easy to de- 
termine the RG flow equations for the magnetic cou- 
pling constants. Our parameter space is spanned by 
H = (g, h, 7s, 7t, <?3, h 3 , £3), and our approximations are 
valid only in the vicinity of the fixed point (FP) fi* = 



(9*,h*,~/*nZ,9th* 3 ,kZ), with l/g* 3 



0, and 



g*, h*, 7*, and 7 t * the FP values of these coupling con- 
stants for the magnetic-field universality class of non- 
magnetic electronsE We therefore immediately linearize 
about this FP. With f3 3 = l/g 3 , and I = 1/k the RG 
length scale, we find 



dl 



dh 3 
dl 

dk 3 
dl 



e- 9 -\r 11 -2{r s +r t \ 

2h* 3 ' 



f22\ 



T /2*3 k 3 



+^ [2(iJ + £)-/£,] fc 3 



(3.21a) 



(3.21b) 



(3.21c) 



where /ft ee /n(7«,7t), ctc - 

The fixed point values that enter Eqs. (3.21) depend 
on whether we consider the long-ranged Coulomb inter- 
action between the electrons, or a short-ranged model 
interaction. We consider here the former, more realistic, 
case. Then we haveo 



5 * = 2e/(l-ln2) 



7* 



-Is 



1 



-2/6 . 

(3.22) 



With this input, we obtain the foll owing eigenvalues for 
the linearized flow equations, Eqs. (3.21), 



Ai = -e/2(l - In 2) + 0(e 2 ) < 
A 2 = -1/(1 -In 2) + 0(e) <0 , 



(3.23a) 
(3.23b) 

(3.23c) 



We see that all three eigenvalues are negative, and the 
fixed point is therefore stable. 



E. Critical behavior 

As we have seen in the previous subsection, the crit- 
ical fixed point for the transition under consideration is 
the same as the one found before for the metal-insulator 
transition of nonmagnetic electrons in the presence of an 
external magnetic field. 3 The asymptotic critical behav- 
ior is therefore also the same. Choosing the correlation 
length exponent v, the critical exponent for the density 
of states P, and the dynamical critical exponent z as the 
three independent exponents, we thus haveo 



v=l/e + 0(l) , 
13= l/2e(l-ln2) 

z = d 



(3.24a) 
(3.24b) 
(3.24c) 



The critical exponent for the conductivity, s = v(d — 2), 



s = 1 + 0(e) 



(3.24d) 



In contrast to the asymptotic critical behavior, the cor- 
rections to scaling are different from any previously stud- 
ied universality class for metal-insulator transitions. The 
reason for this is the presence of the three irrelevant op- 
erators 1/(73, h 3 , and k 3 in our model. We will not go 
through a complete analysis of the corrections to scaling 
here, but only mention that they lead to a nonanalyticity 
in the magnetization as one crosses the metal-insulator 
transition, even though the magnetization is not critical. 
To see this, we recall that the magnetization is propor- 
tional to a frequency integral over the inverse of the one- 



point vertex T^, see Eq. ( 3.5b ) above and Eq. (2.7c) in 
(I). The extra frequency integration makes the integral 
finite for all d > 0, and the one-loop contribution to the 
magnetization is simply proportional to I/53, h 3 , and k 3 . 
Since A3 has the smallest absolute value of the three neg- 
ative eigenvalues given in Eq. ( p3| ), the magnetization 
at T = behaves like 



M(t,T= 0) cx const. + t 



-1/A3 



(3.25a) 



where t is the dimensionless distance from the critical 
point. At criticality as a function of temperature we have 



M (t = 0, T) cx const. + T~ Xa/z 



(3.25b) 



Putting e = 1 in our one-loop approximation yields 
—\ 3 /z = 0.086.... Our theory thus predicts that the 
metal-insulator transition is reflected in the magnetiza- 
tion in the form of a very slow temperature dependence. 

More generally, the existence of very slow corrections 
to scaling indicates that it will be very difficult, if not im- 
possible, to observe the true asymptotic critical behavior 
at the ferromagnetic metal-insulator transition. If we ex- 
trapolate our one- loop results to d = 3 by putting e = 1 , 
we have 1^3 —0.26. This mean that in order to obtain 
the critical exponents with an accuracy of 10% one needs 
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to be within about 0.01% of the critical point, t < 10~ 4 . 
This is not achievable for any metal-insulator transition 
observed so far.0 Any observed critical behavior at larger 
values of t will yield effective exponents that contain con- 
tributions from the dominant irrelevant scaling variables. 
We note in passing that the same conclusion holds for 
the Anderson-Mott transition_of paramagnetic electrons 
in an external magnetic field.u 



IV. DISCUSSION 

Our chief result is the prediction that the metal- 
insulator transition from a ferromagnetic metal to a fer- 
romagnetic insulator is in the same universality class as 
the one from a paramagnetic metal to a paramagnetic 
insulator in the presence of an external magnetic field. 
It is important to note that this statement holds inde- 
pendent of what the actual critical exponents, which we 
can determine only to lowest order in a 2 + e expansion, 
are in three-dimensional systems. It is also independent 
of the fact that we have considered explicitly only the 
perturbative stability of the nonmagnetic fixed point. In 
principle, the full flow e quatio ns that follow from the one- 
loop calculation in Sec. Ill B could contain other critical 



fixed points. This is a question that remains to be inves- 
tigated; here we just mention a possible scenario. 

In both the magnetic field and ferromagnetic material 
cases a different universality class for the MIT is easy to 
envisage. First note t hat t he existence of two conduc- 
tivities, ct ± (cf. Eqs. (2.16)), or equivalently, two diffu- 
sivities, just reflects the fact that either a magnetic field 
or a finite spontaneous magnetization leads to a split- 
ting of the energy band. The subband with fewer (more) 
electrons, that have spins aligned in (opposite to) the 
direction of the magnetic field or spontaneous magneti- 
zation is called the minority (majority) subband. If the 
magnetic energy scale is large compared to other interac- 
tion energy scales and comparable to the Fermi energy, 
then the two subbands are well separated. A polariza- 
tion scenario for the MIT is that the minority subband 
carriers become localized first and then act as_a static 
random field for the majority mobile carriers Ej In this 
scenario, the MIT occurs when the carriers in the major- 
ity band become localized. In this case the MIT is one for 
spin-polarized, or effectively spinless, electrons. This is 
mathematically described by the so-called singlet-only or 
magnetic impurity universality class that was discussed 
One might thus expect a multicritical point 



in Ref. 21 



separating the magnetic field universality class, which 
was discussed above and is relevant for small values of 
the magnetization or the magnetic field, from the singlet- 
only or polarization universality class at large values of 
the magnetization. Experimentally, such a multicritical 
point could be probed by increasing the magnetic field in 
the case of an MIT in an external magnetic field, or by 
effectively increasing the triplet interaction for a spon- 



taneously magnetized system. Theoretically, it remains 
to be seen whether such a behavior is described by our 
complete flow equations. This point will be investigated 
in a future publication. 

In any event, it would be very interesting to compare 
experiments on a ferromagnetic metal-to- insulator transi- 
tion, which has not been studied so far, with the existing 
results for nonmagnetic systems in a magnetic field.EI The 
equivalence of the two universality classes also leads to 
the conclusion that the existing theory for the nonmag- 
netic transition in a magnetic field is incomplete since it 
misses important corrections to scaling. 

From a theoretical point of view, this result is a priori 
rather surprising. The critical behavior at the metal- 
insulator transition, and hence the universality class, is 
determined by the structure of the soft modes in the sys- 
tem, at least near two-dimensions. Since ferromagnetism 
leads to additional soft modes, namely the Goldstone 
modes or spin waves, compared to paramagnetic metals, 
one would expect the critical behavior to change. The 
reason why it does not lies in the fact that the Goldstone 
modes do not lead to a singular correction to the con- 
ductivity in d = 2, in contrast to the diffusive soft modes 
that are also present in the absence of ferromagnetic long- 
range order. Since these singular corrections drive the 
transition in low dimensions, and since the Goldstone 
modes are the only substantial difference between the 
soft-mode spectra of ferromagnetic systems and param- 
agnetic systems in a magnetic field, respectively, the fact 
that the universality class remains unchanged is at least 
plausible. 

We finally mention again that one of our motiva- 
tions for the present study had been the observed ap- 
parent metaljinsulator transition in certain 2 — d elec- 
tron systemsja which contradicts the results of orthodox 
theories and is not understood. Since it is known that 
ferromagnetic fluctuations enhance the conductivity in 
d = 2,0 it was a plausible hypothesis that ferromagnetic 
long-range order might have an even stronger effect and 
lead to a metallic phase in d = 2. Our results rule out 
this possibility, at least on a perturbative level. 
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